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Lovesey and Staub have argued [S.W. Lovesey and U. Staub, Phys. Rev. B 61, 9130 (2000)] that 
experimental data on the temperature dependence of the paramagnetic relaxation of lanthanide 
ions doped into YBa2Cu306+a: are in agreement with the predictions of a model that describes the 
relaxation as due to the scattering of phonons via a magnetoelastic interaction. By generalising 
their model I show that the level of agreement is strongly dependent on the number of intermediate 
lanthanide energy levels included in the calculation, and that inclusion of a more complete set of 
levels leads to very different results that do not necessarily support the phonon damping picture. 



PACS numbers: 74.72.Bk, 74.25. Ha, 76.20.+q, 76 

In a recent papeiE, Lovesey and Staub considered how 
lattice vibrations could contribute to the relaxation of a 
transition between the crystal field split energy levels of 
a paramagnetic ion. Their-work was prompted by a num- 
ber of recent publicationsou reporting neutron spectro- 
scopic measurements of the temperature dependence of 
crystal field transition linewidths for various lanthanide 
ions doped into cuprate superconductors. Lovesey and 
Staub's contention is that the relaxation observed for 
these materials is dominated by decay of the crystal field 
excitation into lattice vibrations, rather than into spin 
excitations on the CuC>2 planes as has been assumed by 
other authors. Lovesey and Staub arrive at this view- 
point by comparing experimental data with the results 
of a model that assumes a magnetoelastic interaction be- 
tween the paramagnetic ion and the lattice. 

The main purpose of this Comment is to raise aware- 
ness of an approximation made by Lovesey and Staub in 
calculating the relaxation from the magnetoelastic cou- 
pling model. The approximation involves the neglect of 
all but 3 of the crystal field split 4/ levels in the re- 
laxation calculation. For the case of Ho 3+ impurities 
in superconducting YBa2Cri3 06+x I will show that the 
use of this 3-level model has led the authors to obtain 
unreasonably good agreement between the calculated re- 
laxation and the experimental data, and that inclusion of 
other thermally populated lanthanide energy levels gives 
results that do not differ significantly from those calcu- 
lated with the exchange coupling model used previously. 
I will also show that the value of the coupling constant in 
the magnetoelastic model derived from the data depends 
strongly on the number levels included in the calculation. 
I will argue, therefore, that the accord between theory 
and experiment presented in Ref. |l| cannot be taken as 
evidence that phonons rather than spin excitations are 
the most important source of relaxation. 

I begin by outlining the factors which enter into a cal- 
culation of the relaxation of crystal field transitions. Fol- 
lowing Lovesey and Staub let us consider a transition 
between two eigenstates \a) and \b) of the lowest energy 
J multiplet of a lanthanide ion in a crystal field. For 
simplicity let \a) be the ground state. The lifetime of 
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the | a) — > | b) transition is finite because of coupling to 
electronic or vibrational degrees of freedom in the host 
system. Relaxation from \b) back to \a) can take place 
either directly or indirectly, the latter proceeding via an 
intermediate lanthanide energy level I7). If I7) is higher 
in energy than \b) then for the indirect transition to take 
place the 4/ ion must be able to absorb energy from the 
host, and this can be done if the temperature T is large 
enough to generate a significant thermal population of ex- 
citations in the host with energies equal to the difference 
in energy between \a) or \b) and I7). The density of host 
excitation states at the relevant energy will clearly be rel- 
evant in determining the importance of a given relaxation 
channel, whether direct or indirect, and a further factor 
is the size of the transition matrix elements connecting 
the states involved. The relevant matrix elements for the 
models under consideration are those of the quadrupolar 
operators Q v in the case of magnetoelastic coupling, and 
of the total angular momentum operators J a in the case 
of exchange coupling. Here v is a symmetry label, and 
a — x,y, z. 

In principle there are as many indirect relaxation chan- 
nels as there are intermediate states I7), but because 
there often exists a wide variation in the sizes of the tran- 
sition matrix elements the effectiveness of each channel 
varies considerably. In Ref. |l] it was argued that inclu- 
sion of only 3 states (\a), \b) and a single I7)) allows the 
relaxation calculation to be simplified while retaining all 
the essential features. 

On the basis of this 3-level approximation Lovesey and 
Staub derived an expression (Eq. (6.7) of Ref. @) for 
the temperature dependence of the transition linewidth, 
and fitted their model to the experimental linewidth 
data for Tb 3+ , Ho 3+ , and Tm 3+ in superconducting 
YBa2Cu 3 06+2; assuming a Debye density of states to de- 
scribe the lattice vibrations. The calculation for Ho 3+ 
is reproduced in_ Fig. 1 together with the data from 
Boothroyd et alB. The energy of the \a) — > \b) transi- 
tion is 0.5 meV, and the value of the intermediate level 
I7) used to generate the calculated curve is 11.8 meV. 

From the quality of the agreement shown in Fig. 1 
it is tempting to be satisfied with the use of the 3-level 
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model, and indeed for the Tb 3+ , and Tm 3+ systems there 
are sound physical arguments based on the energies of the 
excited crystal field levels and size of the quadrupole ma- 
trix elements to suggest that the totality of intermediate 
levels can reasonably be approximated by a single level. 
These arguments are discussed in Ref. In the case 
of Ho 3+ , however, the validity of the 3- level approxima- 
tion is considerably harder to justify because there are 
no fewer than 6 crystal field levels in the energy range 1- 
12 meV several of which have substantial quadrupole ma- 
trix elements connecting them to levels \a) or \b). Hence, 
these levels will contribute to the measured relaxation of 
the 0.5 meV \a) — > \b) transition through indirect pro- 
cesses. 

These considerations make it worthwhile extending the 
magnetoelastic model to include an arbitrary number of 
intermediate levels. For simplicity we will take level |o) 
to be the ground state and level \b) the first excited state, 
and re-label these states as |0) and |1) respectively. In 
this simplification all other intermediate levels I7) are 
higher in energy than |1), as occurs in the experimental 
systems considered here, but it is straightforward to treat 
other sequences of levels. Let ui 1 be the energy of I7) 
relative to the ground state (h — 1). Generalising the 
results of Ref. [l] we can express the linewidth (half width 
at half maximum) in the magnetoelastic model as 
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where, 
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Z v {ui) is the density of phonon modes of symmetry v at 
energy oj , n(u>) = l/(e 0LO - 1), and [3 = l/k B T. c u , c' v 
and c" are proportional to the magnetoelastic coupling 
constants and depend upon the symmetry of the modes. 

The corresponding expression for the linewidth due to 
damping by spin excitations mediated via an isotropic 
exchange interaction Jg x is 
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Here, x"( w ) is the imaginary part of the local dynamical 
susceptibility at the position of the paramagnetic ion. To 
model the normal state of the superconductors we take 
x"(w) oc cu for reasons discussed in Ref. |[ 

With Eqs. (0) and (§ written in this form it is easy to 
see that there exists a simple mapping between Tme and 
Tex, as was pointed out in Ref. [j]. Eq. (0) becomes Eq. 
(||) if we replace the c 2 constants by 2 J| x , the Q v by the 
J a operators, and the function Z v (uS)/oj by "x!'{ui). 

Having stated the basic formulae I will now proceed 
to a comparison of the temperature dependence of the 
linewidth predicted by the two models for the case of 
Ho 3+ in YBa2Cu3 06+x- For this purpose I will make the 
same simplifications as made in Ref. [H, namely to use the 
Debye density of states Z(u>) = 3w V w r> with = 80 
meV for Z v (ui), independent of and to treat all the 
magnetoelastic coupling constants as equal and mode- 
independent. The quadrupole matrix elements are calcu- 
lated from the Ho 3+ wavefunctionSj-derived from the es- 
tablished model for the crystal field!!]. This means there 
is only one unknown parameter c which can be adjusted 
to match the experimental linewidth data in the tem- 
perature range over which the superconductor is in the 
normal state (T > 100 K). 

Fig. 2(a) displays the curve obtained from the magne- 
toelastic model including the first 8 levels of the crystal 
field split J = 8 ground state of Ho 3+ together with the 
same data as on Fig. 1. In principle one could include 
all 2 J + 1 = 17 levels in the calculation, but we prefer to 
apply a cut-off after level 8 (11.8 meV) because there is a 
large gap to level 9 (at 58 meV) ai. which energy the mea- 
sured phonon density of statesEil is much smaller than 
that calculated from the Debye model. This, together 
with the small thermal occupancy at energies of 58 meV 
and above, means that the contribution to the relaxation 
from levels 9-17 will in reality be negligible. A com- 
parison of the calculated curves in Fig. 1 and Fig. 2(a) 
reveals a very significant difference between the 3-level 
and the 8- level approximations. The shape of Tme(T) 
calculated in the 3-level approximation has too much cur- 
vature in the temperature range 10-100 K relative to the 
more complete calculation. The reason for this differ- 
ence, as anticipated earlier, is the neglect in the 3-level 
approximation of 5 of the 6 levels in the energy range 1— 
12meV. These levels become thermally accessible as the 
temperature is raised and cause a significant amount of 
relaxation through indirect transitions. 

What is clear from Fig. 2(a) is that the experimental 
points deviate systematically from the calculated curve 
for temperatures below T « 100 K. This deviation is em- 
phasised in Fig. 2(b) by plotting T/Tme, the ratio of the 
experimental linewidths to the theoretical linewidths de- 
rived from the 8- level magnetoelastic model. Figs. 2(a) 
and (b) can be compared with the corresponding results 
for the spin exchange coupling model shown in Figs. 2(a) 
and (d) of Ref. |5| which exhibit a similar deviation. 

The important point to emphasise, therefore, is that 
if only 3 levels are included in the relaxation calculation 
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then the predicted curve follows the experimental points 
very well, as illustrated by Fig. 1, but as more inter- 
mediate levels are included the experimental linewidths 
fall significantly below the calculated curves both for the 
magnetoelastic model and for the spin exchange cou- 
pling model (the deviation being somewhat larger for 
the latter model). The apparent agreement between ex- 
periment and theory suggested by Fig. 1 cannot, there- 
fore, be used as evidence to favour the magnetoelastic 
model over the exchange coupling model as claimedlij by 
Lovesey and Staub. Indeed, if anything, the failure of 
the magnetoelastic model to describe the data depicted 
in Figs. 2 provides further support for the spin exchange 
model because in this scenario the anomalous reduction 
in linewidth below T c is explained by the opening of a 
gap (the superconducting and/or pseudogap) in the elec- 
tronic excitation spectrum of the Cu02 layers. 

The second point I wish to make concerns the mag- 
nitude of the magnetoelastic coupling constant. As in- 
dicated on Figs. 1 and 2(a) the value of c 2 needed to 
fit the experimental data changes from 1.93 meV 3 to 
0.268 me V 3 as the number of intermediate levels is in- 
creased from 3 to 8. This sensitivity to the number 
of levels included in the model is significant if the cou- 
pling constants for different lanthanides are to be com- 
pared. For example, in the magnetoelastic model pro- 
posed in Ref. ^ c is proportional to the Stevens factor a 
of the lanthanide ion, and so the strength of the para- 
magnetic relaxation for one lanthanide ion can be scaled 
onto another. In Ref. |l| linewidth data for Tb 3+ and 
Ho 3+ in superconducting YF^CusOe+z were found to 
be consistent with the predictions of this scaling hypoth- 
esis when analysed with the 3-level model. Mindful of 
the sensitivity of c 2 to the number of intermediate lev- 
els included in the, relaxation calculation I re-analysed 
the linewidth dataQ for Tb 3+ with a 13-level model and 
found c 2 = 11.5 me V 3 (the equivalent value from the 
3-level model is 24meV 3 ). My estimate for the ratio 
c 2 (Tb)/c 2 (Ho) is 11.5/0.268 = 43, which compares with 
the scalhig-prediction a 2 (Tb)/a 2 (Ho) = 20.7, a factor 2 
differenceEj. Hence, the validity of the claim in Ref. |l| 
that the scaling of the linewidths supports the magne- 
toelastic model for the paramagnetic relaxation is once 
again seen to be dependent on how many crystal field 
levels are included in the calculation. 

I will finish with some brief comments on the two mech- 
anisms, magnetoelastic or exchange coupling, proposed 
to explain the paramagnetic relaxation of lanthanides in 
YBa2Cu30g+a;. It has not been the intention here to crit- 
icize the magnetoelastic model itself. Indeed, the theory 
provided by Lovesey and Staub, summarised in Eq. (|l|) 
above, is an important contribution that permits progess 
to be mad£_in the disentanglement of different sources of 
relaxatiorO. Rather, the purpose of making this Com- 
ment is to illustrate the pitfalls of comparing a minimal 
model with experimental data. The inclusion of more re- 
laxation channels has dramatically changed the quality 
of agreement with the data, and one cannot accept the 



analysis in Ref. [lj as evidence in favour of phonon damp- 
ing. My opinion is that a true assessment of whether 
phononic processes provide an important mechanism can 
only be made when more is known about the phonon 
density of states and magnetoelastic coupling strengths 
for individual phonon modes. 

In my opinion one of the strongest arguments for 
the exchange coupling interpretation is the evidence in 
many of the measurementsEr&t] for an anomalous re- 
duction in linewidth at .a, temperature at or above T c , 
an effect first observedta in the non-cuprate system 
Lai-^Tb^A^. Such an anomaly cannot be produced by 
phonon dampmg. The absence of any anomaly in other 
measurementsQw can perhaps be explained by the size of 
the effect. In most cases where the crystal field transi- 
tion probed is below 1 meV the reduction in linewidth is 
found to be typically 0.02 meV, e.g. Fig. 2(a). This is 
smaller than the error bar on the data around T c in Refs. 

and ^| which are for quasielastic transitions, and so it 
would not have been possible to observe such an anomaly 
with the experimental sensitivity available. 

It is hoped that the continuing debate over phonon ver- 
sus spin fluctuation damping will stimulate a better un- 
derstanding of paramagnetic relaxation mechanisms for 
lanthanide ions in correlated electron systems, and that 
the formulae summarised here may facilitate progress in 
this area. 
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FIG. 1. Temperature dependence of the intrinsic linewidth 
of the ground-state-to-first-excited-state (0.5 meV) crystal 
field transition of Ho 3+ in Hoo.i Yo.9Ba2Cu307. The points 
are the experimental data from Ref. |E] and the line is calcu-i 
lated from the magnetoelastic model of Lovesey and StaubtJ 
with inclusion of 3 crystal field levels. 
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FIG. 2. (a) The same plot as shown in Fig. 1 but with 8 
crystal field levels included in the magnetoelastic model cal- 
culation and a reduction in the c 2 coefficient from 1.93 me V 3 
to 0.268 meV 3 . (b) Temperature dependence of the reduced 
linewidths obtained by dividing the data by the theoretical 
curve in (a). 
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